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HEAVY-TRAFFIC LIMITS FOR WAITING TIMES IN 
MANY-SERVER QUEUES WITH ABANDONMENT 1 

By Rishi Talreja and Ward Whitt 

Columbia University 

We establish heavy-traffic stochastic-process limits for waiting 
times in many-server queues with customer abandonment. If the sys- 
tem is asymptotically critically loaded, as in the quality- an d-efficiency- 
driven (QED) regime, then a bounding argument shows that the 
abandonment does not affect waiting-time processes. If instead the 
system is overloaded, as in the efficiency-driven (ED) regime, fol- 
lowing Mandelbaum et al. [Proceedings of the Thirty-Seventh An- 
nual Allerton Conference on Communication, Control and Comput- 
ing (1999) 1095-1104], we treat customer abandonment by studying 
the limiting behavior of the queueing models with arrivals turned 
off at some time t. Then, the waiting time of an infinitely patient 
customer arriving at time t is the additional time it takes for the 
queue to empty. To prove stochastic-process limits for virtual wait- 
ing times, we establish a two-parameter version of Puhalskii's invari- 
ance principle for first passage times. That, in turn, involves proving 
that two-parameter versions of the composition and inverse mappings 
appropriately preserve convergence. 

1. Introduction. In the past, heavy-traffic stochastic-process limits for 
waiting times in queueing systems have been proven using Puhalskii's invari- 
ance principle for first-passage times [11] together with established stochastic- 
process limits for queue-length processes; e.g., [1, 4, 13, 19]. For instance, for 
an n server system, denoting the arrival, departure and number-in-system 
processes by A, D and X, respectively, the virtual waiting time at time t 
can be represented as the following first-passage time: 

(1.1) V(t) = inf{s >0\D(t + s)> X(0) + A(t) - (n - 1)}. 
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When one attempts to incorporate an abandonment process L in the model, 
one cannot simply write the virtual waiting time at time t (the waiting time 
of an infinitely patient customer arriving at time t) as 

(1.2) V(t) = inf{s >0\D(t + s) + L(t + s) > X(0) + A(t) - (n - 1)}, 

because the term L(t + s) may include customers that arrived to the system 
after time t and then abandoned. 

When the system is asymptotically critically loaded in fluid scale, as in 
the quality-and-efficiency (QED) regime, we show that abandonments do 
not affect fluid and diffusion limits for the virtual waiting time. We do this 
by bounding the virtual waiting time process from above and below by 
processes of the form (1.1) and (1.2), respectively. These bounds are shown 
to be asymptotically equivalent, giving us a limit for our virtual waiting 
time process. For the Markovian M/M/n + M model in the QED regime, 
the virtual waiting time limit was first stated in Theorem 3 of Garnett et al. 
[3], but the proof given there is not correct. We give a correct proof here. 

When the system is not asymptotically critically loaded in fluid scale, we 
must be more careful. One way to address this problem, which was first 
suggested by Mandelbaum et al. [7], is to consider the system with the 
arrival process stopped at a fixed time t, so that the abandonment term 
in the first-passage time expression (1.2) does not include customers that 
arrive to the system after time t. We can then use the invariance princi- 
ple for first-passage times to get the desired limit for the fixed t. To prove 
the full stochastic-process limit for virtual waiting times, however, we need 
a corresponding two-parameter limit for the queue-length process with ar- 
rivals turned off and a two-parameter version of the invariance principle for 
first-passage times. We establish these results here. We establish the two- 
parameter version of Puhalskii's invariance principle for first-passage times 
by proving that two-parameter versions of the composition and inverse map- 
pings are convergence-preserving. This portion of our work is along the lines 
of [14] and Chapter 13 of [15]. 

We apply the two-parameter version of Puhalskii's theorem to establish 
many-server heavy-traffic stochastic-process limits for waiting times in the 
M/M/s + M model in the efficiency-driven (ED) regime. As indicated above, 
the line of reasoning here was initiated by Mandelbaum et al. [7], which 
builds on the strong-approximation approach in [6] . In [7] the need to turn off 
the arrival process to properly treat customer abandonment was recognized, 
but the two-parameter version of Puhalskii's theorem needed to complete 
the argument was not stated or proved. The statement and partial proof 
were developed in an unpublished manuscript [8]. The most difficult step in 
proving the ED result is not establishing the two-parameter generalization of 
Puhalskii's theorem, but is instead establishing the stochastic-process limit 
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for the queue-length process with the arrival process turned off, Proposition 
6.2 here. Mandelbaum et al. [8] establish corresponding results in the more 
general setting of "Markovian service networks" [6], but that part of their 
work remains to be completed. Our task is made easier by considering only 
the M/M/s + M model. Nevertheless, we believe that the results and line 
of reasoning here apply to more general models. 

We also establish the heavy-traffic limit for the steady-state waiting time 
in the ED regime by a different argument. As should be anticipated, it is a 
normal distribution centered at the fluid limit. We show that the variance 
has a simple form, which is unclear a priori. Our result is a special case 
of the ED limit for the MjMj s + GI steady-state waiting time in (6.12) 
of Zeltyn and Mandelbaum [20] . Our proof is interesting because we use a 
stochastic-process limit to prove the one-dimensional limit. For the QED 
regime, the limit for the steady-state waiting time distribution is given in 
Theorem 3 of [3]. 

Notation. For a,6,c£l, let a\/b = max{a,b}, aAb = min{a,&}, [c] b a = 
max{min{c, b}, a}, and a + = a V 0. 

For a complete separable metric space (CSMS) (S,d), we let C([0, oo),<S) 
denote the set of continuous functions from [0, oo) to S and define C = 
C([0, oo),R) and Cq = C([0, oo),C). We endow these spaces with the topol- 
ogy of uniform convergence over compact sets (u.o.c). For x S C([0,oo),5), 
let || a? — y\\t = sup 0<s<t d(x(s) , y(s)) . A metric m inducing the topology of 
u.o.c. is 



Notice that m(x n , x) — > as n — > oo if and only if ||x n — x\\t — > as n — > oo 
for all t > 0. 

For a CSMS S, let D([0, oo),<S) denote the set of right-continuous func- 
tions from [0, oo) to S having left limits everywhere on (0, oo) (see [2] and 
[15]) and let D = D([Q, oo),M) and D D = D([0, oo), D). Let e denote the 
identity element in D, i.e., e(t) = t for t > 0, and let r\ denote the constant 
function 1, i.e., rj(t) = 1 for t > 0. We use (D,T) to denote the space D 
with the topology T. We use 71, 7^) to denote the space Dp with 
the topology T\ on the outside D space and the topology T2 on the inside 
D space. We use J\ and M\ to denote Skohorod's J\ and M\ topologies, 
respectively, and U to denote the topology of u.o.c. as defined above. Note 
the spaces (Dp, Ji,M\) and (Dp, Ji, Ji) are well-defined topological spaces, 
but (Dp, Mi, Mi) and (Dd,M\, Ji) are not, since the M\ topology has only 
been defined on the spaces Z)([0,oo),lR fc ), for k > 1 (see Chapter 12 of [15]). 



(1.3) 
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We can identify functions in Cc and Dp as two-parameter functions x = 
x(-,-). For x G Dp and s > 0, we will often find it convenient to refer to 
x(s, •) G -D([0, oo),R) as x s . We also work with product spaces such as D 2 D = 
Dp x Dp, where we assume the appropriate product topology. For x G D, let 
= sup < s <5 |x(s)| and for x G Z?r>, let ||x||s,t = sup < s < 50 < t < T \x(s,t)\. 
(This is consistent with our notation \\x — y\\t used in (1.3); since S = K, 
we now have a norm.) We will mainly be considering limit processes with 
continuous sample paths. It is thus significant that for x G Cc convergence 
x n — ► x in (Do, Ji, Ji) or (Z?d, J±,Mi) is equivalent to — ► x in (Dp, U, U), 
which, in turn, is equivalent to \\x n — x\\st ~^ f° r an S,T > 0; see Theorem 
2.1 below. 

Let => denote convergence in distribution. For a sequence of M-valued 
random variables (X n ) n >i, we write X n = op(/(n)) if 

- " =4> in 1 as n->oo. 
/(») 

Organization. In Section 2, we begin by proving various functions in 
Dp are convergence preserving. Some of these results are given in greater 
generality than needed here. The additional generality may be useful in 
future work. In Section 3, we use a simple bounding argument to prove 
heavy-traffic limits for virtual waiting time processes when the abandonment 
process is asymptotically negligible in fluid-scale. In Section 4, we use the 
results in Section 2 to prove heavy-traffic limits for virtual waiting time 
processes assuming associated limits hold for basic queueing process with the 
arrival process stopped. In Sections 5 and 6 we use the results in Sections 3 
and 4, respectively, to prove waiting-time stochastic-process limits for the 
M/M/n + M model in the QED and ED limiting regimes. In Section 6.2, we 
establish the heavy-traffic limit for the steady-state virtual waiting times in 
the ED regime. Section 7 contains the proofs of lemmas used in Sections 3 
and 6. 

2. Functions in Do. In this section, we prove continuity results for func- 
tions in Dp. We focus on composition and inverse in Sections 2.1 and 2.2, 
respectively. We use the results in Sections 2.1 and 2.2 to extend Puhalskii's 
invariance principle for first-passage times [11] to functions of two parame- 
ters in Section 2.3. We prove continuity results for integral mappings in Sec- 
tion 2.4 and for a projection mapping in Section 2.5. For our two-parameter 
results, the space can be either (Dp, Ji, J\) or (Dp, J±, Mi), unless the Mi 
topology is explicitly assumed on the inside space, but this usually does not 
matter because our limits will usually be assumed to be in Cc- 

Let Z?2,o he the subset of functions x in Dp such that x(0, 0) > 0. Let 
and L>2, n be the subsets of functions in -D2,o that are nondecreasing and 
strictly increasing in each coordinate, respectively. Similarly, define C2,o = 
Cc H -t>2,(b C2,t = Cc H D>2 : -\ and C^-fj = Cc H ^2,TT- 
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2.1. Composition. First, we show that, as in the one-parameter case, 
convergence in to a point in Cc is equivalent to uniform convergence 
over compact sets. We need the following lemma, whose proof is an easy 
exercise in analysis. 

Lemma 2.1. Let (S,d) be a metric space. Let (x n ) n >i and x be functions 
from a compact set EcR 4 to S , k > 1, and let x be continuous. Then 
sup tgE d(x n (t),x(t)) — > as n — > oo if and only if for all sequences (i n )n>i Q 
E such that t n —*t as n— > oo we have x n {t n ) — > x(t) in S as oo. 

Theorem 2.1. LfxeCc, thenx n —>x in (D D , Ji, J{) or (D D , Ji,M x ) 
as oo if and only if x n — > x in (Dd, U,U) as n — > oo. 

Proof. Let T be J\ or M\. Since x G Cc the convergence x n — > x in 
(De>, J\,T) is equivalent to convergence in (Dd,U,T). Therefore, for all 
T > 0, we have 

(2.1) sup dfix^^x 1 ) — > asn^oo, 

te[0,T] 

where dq- is a metric on D inducing the topology T. By Lemma 2.1, (2.1) 
holds if and only if for every T > 0, t G [0, T] , and sequence (i n )n>i ^ [0>^1 
such that t n — > t as n — > oo we have — > x* in (D([0,T],]R),T). But since 
i 1 £ C for each t e [0,T] this is equivalent to x^ x t in (D([0,T],M), U) 
for all T > 0. Again, applying Lemma 2.1, this convergence holds if and only 
if for for each S, T > 0, (t, s) G [0,T] x [0,5], and sequence ((t n , s n )) n >i C 
[0, T] x [0, S] such that (t n ,s n ) — > (t, s) as n — > oo we have x n (t n ,s n ) — > x(t, s) 
in R as n — > oo. But then again by Lemma 2.1, this is true if and only if 
x n — * x in (Dd, U, U). □ 

We now show that addition and multiplication on D 2 D are measurable and 
continuous at points in C^. In both cases, to prove measurability, we use 
Lemma 2.7 of [14], which we state here for completeness. 

Lemma 2.2 (Lemma 2.7 of [14]). The Borel a-field on Dd coincides 
with the Kolmogorov a-field, i.e., the a-field generated by the coordinate 
projections. 

Let addition on D D be defined pointwise: For x,y G Dd, 
{x + y){s,t) =x(s,t) + y(s,t). 

Theorem 2.2 (Continuity of addition in Dd)- Addition on D D is mea- 
surable and continuous at points in Cq. 
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Proof. Measurability follows by Theorem 4.1 of [14]. The argument 
there is based on Lemma 2.2. Continuity follows from Theorem 2.1 and the 
triangle inequality. For each S, T > 0, 

\\(x n + y n ) - (x + y)\\ s ,T < \\x n - x\\s,t + \\y n -y\\s,T- □ 

Similarly, define multiplication on by 

(xy)(s,t)=x(s,t)y(s,t) 

for x, y G Dv. 

Theorem 2.3 (Continuity of multiplication in Do). Multiplication on 
D 2 D is measurable and continuous at points in C c . 

Proof. Again, measurability follows using the same argument as in 
Theorem 4.1 of [14]. Continuity at continuous limits follows from Theorem 
2.1 and the triangle inequality: 

\\x n y„ - xy\\s,T < \\x n y n ~ x n y\\ s ,T + \\x n y ~ xy\\ s ,T 

< \\x n \\s,T\\yn ~ V\\s,T + \\y\\s,T\\ x n ~ x \\s,T- 

□ 

Since x n — > x, sup n>1 ||x n ||s,T < °o, and we have our result. 

Define the composition mapping o 2 from Do x Dp to -D2,t by 
(2.2) {xo 2 y)( s ,t)=x(s,y( S ,t)). 

Notice that y only enters in the second argument of x. As usual, we re- 
quire that the internal function be in -D2,t- This is needed to ensure that 
the composition x o 2 y is an element of Dp (see Example 13.2.1 of [15]). 
Alternatively, we can define the composition mapping as a mapping from 
Cc x Cc to Cc- For each s > 0, we have (x o 2 y)(s, •) = x(s, y(s, ■)) = x s oy s , 
where o is the standard composition map on D. 

Then we have the following extension of Theorem 3.1 of [14] and Theo- 
rem 13.2.2 of [15]. 

Theorem 2.4 (Convergence preservation of o 2 ). 

1. The composition function o 2 mapping Dp x D 2j i- into Do defined in (2.2) 
is measurable and continuous at each (x,y) S Cc x C 2 i-. 

2. The composition function o 2 mapping Cc x Cc into Cc is continuous. 

Proof. For part (1), measurability follows from Lemma 2.2 and the 
fact that the composition mapping from D x to D is measurable (see 
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Theorem 13.2.1 of [15] and page 232 of [2]). For both parts (1) and (2), we 
use the triangle inequality to prove continuity. For each S, T > 0, we have 

\\x n °2Un-X o 2 y\\s,T < \\x n °2 Dn - x o 2 y n \\s,T + \\x o 2 y n - x o 2 y\\s,T- 

Let T' = sup n ||y n ||s,T- Then the first term converges to zero since it is 
bounded by ||x n — »||st' — > 0. The second term converges to zero since 
x G Cc- The mapping in part (2) is measurable because it is continuous. 

□ 

We now give the main result of this subsection. It is an extension of 
Theorem 13.3.3 of [15]. 

Theorem 2.5 (Convergence preservation of o 2 with nonlinear centering). 
Let (x n ,y n ) G (Dd x D 2 ^) U (Cc x Cc) and y G Cc- Let x have continuous 

partial derivative x' with respect to its second parameter, i.e. x'(s,t) = dX g t '^ , 
and let c n — > oo . If 

c n (x n -x,y n -y)^(u,v) in D 2 D , 
where u,v G Cc, then 

c n (x n °2 Vn ~ x o 2 y) -> u o 2 y + (x' o 2 y)v in D D . 

PROOF. We follow the proof of Theorem 13.3.3 of [15]. Let (x n ,y n ) G 
Dd x £*2,t; ^ e °th er case can be proven analogously. Note that 

C n (x n o 2 y n -X o 2 y) = c n (x n - x) o 2 y n + c n (x o 2 y n -X o 2 y). 
By our assumptions, we have 

[cn{x n ~ x),c n (y n -y),y n ) -* (u,v,y) in D%, 
so, applying Theorems 2.2 and 2.4, gives us 

(2.3) (c n (x n o 2 y n -x o 2 y n ) + (x' o 2 y)c n (y n -y))^uo 2 y+ (x o 2 y)v. 
Now, for all S, T > 0, 

\\c n (x n o 2 y n -x o 2 y) - c n (x n o 2 y n -x o 2 y n ) - c n {x! o 2 y)(y n - y)\\ S) T 
< \\c n {x °2Vn-x o 2 y) - c n {x' o 2 y)(y n - y)\\ ST 

(2.4) 

rvn(s,t) 

■ sup 

0<s<5 
0<*<T 



x'(s, u) du - c n x'(s, y(s, t))(y n {s, t) - y(s, t)) 
y(s,t) 



< sup 

0<s<5 
0<t<T 



sup x'(s,u)-x'(s,y(s,t)))c n (y n (s,t)-y(s,t)) 

y(s,t)<u<y n (s,t) ' 



The last expression goes to zero since x' is uniformly continuous over bounded 
intervals in its first parameter and c n (y n — y) — ► v in Dd. Now, combining 
(2.3) and (2.4) and using the triangle inequality gives us our result. □ 
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2.2. Inverse with centering. We now move on to prove results about 
inverse maps with linear and nonlinear centering in Dp. The inverse we 
consider here is the usual inverse applied to only the second argument. Let 
Z?2,u be the subset of functions x in Dd such that x(0, 0) = and x(s, •) 
is unbounded above for each s > 0. As before, let D2 jU ,i = D2,u H D% t and 
Z) 2 ,«,TT = ^2,ii H Z?2,tT • We define the inverse map on the subset Z?2,« of D 
as follows. For x G -D2,u> let the inverse of x be 

x~ 1 (s, i) = inf{u > | x(s, u) > t}. 

Notice that for each s > 0, 

(x' 1 ) 3 = x~ 1 (s, •) = inf{u | x(s,u) > ■} = inf{u \ x s (u) > •} = 

where the inverse on the right is the standard inverse on D. 

We then have the following result, which is analogous to Theorem 13.6.1 
of [15], which is proved in Section 7.6.1 of [16]. The theorem exploits the M\ 
topology on the inner space D. Note the second condition of the theorem 
holds if x £ Cc, which will be the case when we apply the theorem below. 

Theorem 2.6 (Convergence preservation of inverse). The inverse map 
on (D2 tU ,Ji,Mi) is measurable and continuous at all x G -t>2,« such that: 

1. x- 1 (s,0)=0/or alls>0, 

2. x is right- continuous at 0, uniformly in s G [0, S] for each S > 0, i.e., for 
all S > and e > 0, there exists S > such that 

(2.5) sup \x s {5) -x s {0)\ <s. 

0<s<S 

Proof. Measurability follows by Lemma 2.2 and the fact that the in- 
verse map on D is measurable, by the argument in the proof of Theorem 
13.6.1 of [15]. 

Suppose x n — > x in (-D2,u, Ji, Mi). Then, for each S > 0, there exists a 
sequence of increasing homeomorphisms on [0, S] such that 

||An — e||sV sup djvfi (^n> £ A?l ^) — ► as n- > oo. 

0<s<S 

It now suffices to show that (x" 1 )** — > (x -1 ) An ( s ) in (D,Mi), uniformly in 
s € [0,5]. To show this, we can use the same continuity argument as in the 
proof of Theorem 13.6.1 of [15]. In particular, we can construct a sequence 
of functions x* in Dp such that for n large enough (x*) -1 (s,0) = for all 
s G [0,S] and ((x*) s ) _1 is close to (a^) -1 in (D,U) uniformly in s G [0,S]. 
In order for that argument to go through uniformly in s G [0,5], we need 
the condition (2.5). The reason we want (x*)~ 1 (s,0) = for each s G [0,5] 
is that if we have this property then each M\ parametric representation of 
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(x*) s serves as a parametric representation of ((x*) 5 )" 1 with the roles of the 
coordinates switched (see Lemma 13.6.8 of [15]). Therefore, for large enough 

n, 

sup ^((s- 1 )',^- 1 )^) 

0<s<S 

= SUp dMAiO' 1 ,^^)' 1 ) 
0<s<S 

< sup d M MOT 1 ,(x Xn{s) r 1 )+ sup d Ml (((x;)- 1 ) s ,(x^)- 1 ) 

0<s<S 0<s<S 

< sup d Ml (((x* n yr\(x Xn(s) r 1 )+ ^p ii (to -1 )' -to -1 II 

0<s<5 0<s<5 

= sup d Ml (( x * n y,x x ^)+ sup iK^r^-Kr^Ho 

0<s<5 0<s<5 

as n — ► oo. □ 

Remark 2.1. Theorem 7. 1 of [14] , which proves continuity of the inverse 
function in (D,Mi), is incorrect. We need the limit process x to have the 
extra property that x _1 (0) = 0. This has been corrected in Section 13.6 of 
[15]. 

We now prove convergence results for inverse with centering as in Section 
13.7 of [15]. Here, we let the function e 2 € Dp be defined as e2(s,t) = t for all 
s, t > 0. Also, for each T > define the maximum jump function Jt '■ D — > R 
by 

Jt(x) = sup {\x(t) — x(t—)\}. 

0<t<T 

Theorem 2.7 (Inverse with linear centering). Suppose that c n (x n — e 2 ) — > x 
in Dd as n — > oo, where x n € -D 2)n , % S Cc, y(0) = 0, and c ra — > 00. T/ien 
Cnix^ 1 — e 2 ) —x in Do as 00. 

Proof. By the triangle inequality, for each S, T > 0, 

\\c n (x~ 1 -e2) + x\\ s ,T < ||c„(^ 1 -x n o 2 x^ 1 ) + 3;|| SiT + ||c n (x n ,02X^ 1 -e2)||5,T- 

Applying the composition map, o 2 (Theorem 2.4) with the assumed conver- 
gence and x" 1 — ► e 2 (Theorem 2.6) gives us convergence of the first term 
above to 0. Furthermore, by Corollary 13.6.2 of [15], we have 

||c„(x n o 2 x" 1 - e 2 )||s,T = c„ sup \\x s n o (x*)" 1 - e|| T 

0<s<5 

< c n sup J {x s n) -i {T) (x s n ) 

0<s<S 

= sup J(x.)-i(T)(Cn(<-e)). 
0<s<S 
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Since c n (x n — e) — > x and x is in Cc, the last expression above converges to 
as n — > oo. This gives us our result. □ 

We now move on to inverse with nonlinear centering. The proofs of our 
next two results, Theorems 2.8 and 2.9, follow the proof of Theorems 13.7.2 
and 13.7.4 of [15] exactly, with results on composition and inverse maps in 
D there replaced by our new results on composition and inverse maps in 
Djj here. 

Theorem 2.8 (Inverse with nonlinear centering). Suppose that c n (x n — 
A) — > u as n — > oo in Dp, where x n G D2 ;U > u G Cc, u(0) = 0, c n — ► oo, and 
A is absolutely continuous with continuous positive derivative A' with respect 
to its second parameter. Then 

t -i \-U -u°2*~ l . n 

c n {x n — A j — ► — — r 2n as n — > oo, 

A o 2 A 1 

where (u/v)(s,t) = u(s,t)/v(s,t), for s,t>0. 

2.3. Two-parameter version of Puhalskii's theorem. Finally, we have our 
two-parameter version of Puhalskii's theorem. This is similar to Theorem 1 
in [11] and Theorem 13.7.4 in [15], which is an extension to the Mi topology. 
Note that the limits x, y, u and v below are all continuous. Again, we omit 
the proof because it is a direct analog of the proof of Theorem 13.7.1 in [15], 
using the results above. 

Theorem 2.9. Suppose that x n £ D2 ;U > Un G ^2,T; c « ~~ * 00 > 

c n (x n -x,y n -y) -> (u,v) in D 2 D , 

where u,v £ Cc, x G C2 jU and y G C2,T- Furthermore, suppose x is absolutely 
continuous with a continuous positive partial derivative x' with respect to its 
second parameter, where x' G Cc- Then 

-i -i v ~ u °2 °2 y 

(2.6) c n (j; n o 2 y n -x o 2 y)^ — — mD D asn^oo. 

x o 2 x o 2 y 

Remark 2.2. Although we require here that x^ is unbounded for each 
s > and n > 1, the proof goes through as long as the expressions on the 
left-hand side of (2.6) exist as elements of Dp. A sufficient conditions for 
this is sup u>Q x s n (u) > sup u>0 y B n (u) for each s > and n> 1. 

Remark 2.3. If we had defined the inverse process for x G D 2)U and 

y g £> 2 ,T by 

(x _1 o 2 y)(s, t) = inf{n > | x(s, u) > y(s, t)}, 
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then x" 1 o 2 y would not necessarily be an element of Do- However, if y s is 
piecewise constant for each s > 0, then x _1 o 2 y is an element of Do- In fact, 
if we make the assumption that y s n is piecewise constant for all s > and 
n > 1, then the theorem holds with this alternative definition of inverse (see 
the remark in [11]). This will be useful in Section 4. 

2.4. Integral mappings. We now prove continuity and measurability re- 
sults for integral mappings in Dp, which are used to prove the two-parameter 
version of our ED queue-length heavy-traffic limit below. We first state a 
continuity result for the basic integral mapping for functions in Do, where 
the integration is with respect to one of the parameters. It is analogous to 
Theorem 11.5.1 of [15], which is for elements of D. 

Theorem 2.10. Let >M be a uniformly continuous function and 

suppose that x n — > x in Do as n — > oo and x G Cc- Then for the integral 
mapping f : Do — > Do defined by 



we have f(x n ) — ► f(x) in Do as n — > oo and f(x) G Cc- 

Proof. First, we show that f(x) G Cc- For each fixed s > 0, f(x)(s,-) 
is clearly in C. Consider a sequence (sn)n>i 

C [0,oo) such that s n — > s as 

n — ► 00. Then, for each T > 0, we have 



The last expression goes to as n — > 00 by bounded convergence and con- 
tinuity of g. Since f(x) G Cc, by Theorem 2.1 it is enough to show that 
f(x n ) — >■ f(x) in (Do, U, U). For each S, T > 0, we have 



Since x G Cc, we have \\x n — x\\s,t — > as n — > 00 for all S, T > 0. Thus, by 
the uniform continuity of g, the right-hand side of (2.7) goes to as n — > 00. 

□ 

We now give continuity results for regulator mappings that are general- 
izations of the regulator mapping of Theorem 4.1 of [9]. 






(2.7) 
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Theorem 2.11 (Continuity of regulator map in D). Consider the inte- 
gral representation 

x(t) = y(t) + / h(x(u),u)du, t>0, 



o 



where h:Mx [0,oo) — > M satisfies the following properties: 

1. h{0,t) = for allt>0. 

2. For each T > 0, there exist c±,C2 > such that, for each x±,X2 € D ant 
homeomorphism X on [0, T] with strictly positive derivative A', we have 



(21 



o 



\h(xi(u),u) — h(x2(X(u)) , X(u))\ du 



< c\ || A — c\\t + C2 / \xi(u) — X2(X(u))\ du 

for all t £ [0,T]. 



o 



This integral representation has a unique solution x, so that it defines a 
function f :D — ► D mapping y into x = f(y). The function f is continuous 
provided its domain and range are both equipped with either: (i) the topology 
U or (ii) the topology J\. Furthermore, ify&C, then f(y) £ C . 

Proof. That / is a well-defined map from D to D follows as in the 
proof of Theorem 4.1 of [9]. To prove continuity in the topology of uniform 
convergence over bounded intervals, notice that the condition (2.8) implies 
that there exists C2 > such that for any x\, X2 £ D we have 

t rt 

\h(xi(u),u) — h(x2(u),u)\du < C2 / \x\(u) — X2(u)\ du 
o J o 

for all t > 0. The proof then follows using GronwalFs inequality as in the 
proof of Theorem 4.1 of [9] 

We now show that / is continuous when both its domain and range are 
equipped with the J\ topology. We mostly follow the proof of Theorem 4.1 
of [9]. Let T > be a continuity point of y £ D and suppose y n — > y in 
D([0,T],M) with the J\ topology. Then there exists a sequence of increas- 
ing homeomorphisms (A n ) n >i of the interval [0,T] such that ||A n — e||t — ► 
and || y n — y o A n ||r as n — > cxo. It suffices to consider homeomorphisms 
(An)n>i are absolutely continuous with respect to Lebesgue measure 
on [0, T] having derivatives A^ satisfying ||A^ — 1||t — > as n— > oo. Then for 
N large enough so that info< u <T X' n (u) > for all n > N we have 



(t)-x(X n (t))\ 
< \y n (t) - y(X n (t))\ + 



•A„(t) 

h(x n (u),u) du — I h(x(u),u)du 



o Jo 
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<\y n (t)-y(X n (t))\ 

rt ry 
+ / h(x n (u),u)du- h(x(X n (u)),X n (u))X' n (u)du 
Jo Jo 

< \y n (t) - y(X n (t))\ + \\X' n - 1\\ T C \h(x(X n (u)), X n (u))\ du 

Jo 

+ / \h(x n (u),u) - h(x(X n (u)),X n (u))\du 
Jo 

< \y n (t) - y(X n (t))\ + ||A'„ - 1|| t ci||x||tT + ci||A n - e\\ T 
+ C2 / \x n (u) - x(X n (u))\du 



for all n > N and < t < T, where the last inequality follows from (2.8). 
Since \\y n — y o A n ||T, ||A^ — 1||t and || X n — e||r converge to as n — > oo, the 
result follows by an application of Gronwall's inequality. 

Finally, the inheritance of continuity is straightforward as in the proof of 
Theorem 4.1 of [9]. □ 

Note that the continuity of f:(D,Ji) — > (D,J\) established in Theo- 
rem 2.11 implies measurability using the Kolmogorov cx-field. We now extend 
Theorem 2.11 to D D . 

Theorem 2.12 (Continuity of regulator map in Dp). Consider the in- 
tegral representation 

x(s,t) = y(s,t) + / h(x(s,u),s,u)du, t>0, 



o 



where h : M x [0, oo) 2 — > R satisfies the following properties: 

1. h(0,s,t) = for all s,t>0. 

2. For each T > 0, there exist c\,C2 > such that for each si,S2 > 0, x±,X2 € 
Dd, and homeomorphism X on [0,T] with strictly positive derivative X' 
we have 

rt 

\h(xi(si,u),si,u) - h(x2(s2,X(u)),S2,X(u))\du 

< ci\s\ - s 2 | + c 2 1| A - e\\r + c 3 / \xi(si,u) - x 2 (s 2 , X(u))\ du 

Jo 

for all t G [0,T]. This integral representation has a unique solution x, so that 
it defines a function f : Do — > Do mapping y into x = f(y). The function f 
is measurable provided its domain and range are equipped with the topology 
(Do , J\ , J\ ) and continuous provided its domain and range are equipped with 
the topology (Do, U, U). Furthermore, if y E Cq, then f(y) S Cq- 



14 



R. TALREJA AND W. WHITT 



Proof. Consider our integral representation for each fixed s > 0. In 
other words, consider the function / s : Dp — > D mapping y € Dp to the 
solution x £ D of 

x(t) = y(s,t) + / h(x(u),s,u) du. 
Jo 

By Theorem 2.11, f s is a well-defined function that is continuous in the 
topology of uniform convergence over bounded intervals. Therefore, our two- 
parameter integral representation has a unique solution since it has a unique 
solution for each fixed s > 0. 

To show that for y £ Dp, x = f(y) is also in Dp, consider a sequence 
(sn)n>i ^ [0, oo) such that s n I s as n — > oo. Since y G D/j, there exists a 
family of homeomorphisms (A n ) n >i on [0, T] such that ||A n — e\\x — ► and 
||y*n _ y« o A n ||T — > as n — > oo. Again, it suffices to consider homeomor- 
phisms (A n ) n >i that are absolutely continuous with respect to Lebesgue 
measure on [0,T] having derivatives A^ satisfying ||A^ — 1\\t — > as n — > oo. 
Then for N large enough so that info< M <T X' n (u) > for all n > N we have 

\x(s n ,t) -x(s,X n (t))\ 

< \y(s n ,t) -y(s,X n (t))\ 



+ 



h(x(s n , u),s n , u) du 



A„(t) 



h(x(s, u),s, u) du 



< \y(Sn,t) -y(s,X n (t))\ 



+ 







h(x(s n ,u),s n ,u) - h(x(s,X n (u)),s,X n (u))X' n (u)du 





< \y(s n ,t) -y(s,X n (t))\ + || A' 

+ /' 

Jo 



\h(x(s,X n (u)),s,X n (u))\ du 



< \y(s n ,t) -y(s,X n (t))\ 

rt 

' h(x(s n ,u),s n ,u) - h(x(s , X n (u)) , s , X n (u))X' n (u) du 

l||t [ T 

Jo 

\h(x(s n ,u),s n , u) - h(x(s, A n (it)), s, A n (u))| du 

< \y{s n ,t) -y(s,Xn(t))\ + \\X' n - l\\ T C2\\x s \\ T T + a\s n - s\ 

+ c 2 ||A n -e||r + C3 / \x(s n ,u) - x(s, X n (u))\du 
Jo 

for n> N and < t < T. The last inequality follows from conditions 1 and 
2. Since \\y Sn —y s o A n ||T, || X' n — 1||t, \s n — s\ and ||A n — e||r converge to as 
n— > oo, the result follows by an application of GronwalFs inequality. That 
left-hand limits exist for x = f(y) follows using a similar argument. 
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Since the Borel cr-field on Dp coincides with the Kolmogorov cr-field gen- 
erated by the coordinate projections, in order to show measurability, it suf- 
fices to show that the map f s defined above is measurable for each s > 0. 
But f s is simply the composition of the projection y^y s and the regulator 
map ff on D defined in Theorem 2.11 with the function h\ : M x [0, oo) — > M 
given by hf(x,u) = h(x,s,u), for x £ K, u > 0. Since both the projection 
y i— > y s and the regulator map ff are measurable, f s is measurable. Thus, / 
is measurable. 

Continuity in the topology of uniform convergence over bounded inter- 
vals follows using the same argument as in Theorem 4.1 of [9], which uses 
Gronwall's inequality. To show that y £ Cc implies x = f(y) £ Cc is also 
straightforward using Gronwall's inequality. □ 

2.5. Projection mapping. Finally, we will be needing the following pro- 
jection lemma. 

Lemma 2.3. Let x n £ Dd, x £ Cc, and suppose x n — > x in Dq as n — > 
oo. Define y n £ D by y n {t) = x n (t,t) for all t>0. Then y n — > y in D as 
n — > oo, where y(t) = x(t, t) for all t>0. 

PROOF. Apply Theorem 2.1 and the fact that 

\\y n -y\\T= sup \x n (t,t) - x(t,t)\ 

0<t<T 

< sup \x n (s, t) — x(s, t)\ = \\x n — x\\t.t — ► 

0<s,t<T 

as n — > oo. □ 

3. Critically loaded. In this section and the next, we establish heavy- 
traffic limits for virtual waiting times for general G/G/n + G queues, assum- 
ing we have corresponding heavy-traffic limits for the basic arrival, queue- 
length, departure, and abandonment processes, and the fluid limit of the 
queue-length process is identically 1. In Section 5, we apply the main result 
of this section with the basic M/M/n + M QED result, Theorem 5.1, to 
obtain a proof of the QED virtual waiting time result for the M/M/n + M 
model, Theorem 5.2. 

There are a number of interesting waiting time performance measures that 
have been studied in the literature. A customer's potential waiting time is 
the amount of time the customer spends waiting in queue if the customer 
were infinitely patient (see [3]). Therefore, if W is the potential waiting 
time of a particular customer and X is his patience (the amount of time 
he is willing to wait in the queue before being served), then the customer's 
actual waiting time is given by WAX. Another waiting time performance 
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measure of interest is customer waiting time conditional on being served. An 
exact steady-state analysis of this performance measure in the Markovian 
case using Laplace transforms appears in [18], but we do not deal with it 
here. Our main focus is to study virtual waiting time processes, which are 
continuous-time process analogs of potential waiting times. A virtual waiting 
time process is a process (V(t),t > 0), where V(t) is the potential waiting 
time of a hypothetical customer arriving to the queue at time t. 

We assume we are given a sequence of G/G/n + G models. In model 
n > 1, let X n = (X n (t),t> 0), A n = (A n (t),t> 0), D n = (D n (t),t> 0), L n = 
(L n (t),t>0) and V n = (V n (t) , t > 0) , where 

X n (t) = number of customers in the system at time t, 

A n (t) = number of customers that have arrived to the system by time t, 

D n (t) = number of customers that have been served by time t, 

L n (t) = number of customers that have abandoned by time t, 

V n (t) = potential waiting time of a hypothetical customer 

arriving at time t. 



t — ™ t/ = y 

n 

X n = c n (X n — X), A n = c n (A n — A), D n = c n (D n — D), 

L n = c n (L n -L), V n = c n (V n -V), 

where X, A,D,L,V G C and c n — > oo and n/c n — > oo as n — > oo. 
We then have the following result. 

Theorem 3.1. Consider a sequence of G/G/n + G models and suppose 

(3.1) (A n ,D n ,L n ,X n )^(A,D,L,X) mD 4 

as n — > oo, where A, D and L are continuous deterministic functions, D 
has continuous positive derivative D' , L = a] for some c > 0, and X = 1. 
Furthermore, suppose 

(3.2) (A n ,D n ,L n ,X n )^(A,D,L,X) in D 4 

as n — > co where A, D, L and X are continuous processes. Then we have 




in D 



as n — > oo jointly with (3.1) and (3.2). 



For each n > 1 , define processes 

n n n 

and 
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Proof. Recalling the discussion about (1.1) and (1.2), first observe that 
for each t > and n > 1, we can bound the virtual waiting time process from 
above and below by first-passage times: 

(3.3) V*(t) < V n (t) < V:(t), 
where 

Vi{t) = inf{s > 0\D n (t + s) + L n (t + s) 

(3.4) 

> X n (t) + D n (t) + L n (t) - (n - 1)}, 
(3.5) V£(t) = ini{s > Q\D n (t + s) > X n {t) + D n (t) - (n - 1)}. 

For all n> 1, define the first-passage-time processes Z l n = (Z l n (t),t > 0) and 
Z'l = (Z%(t),t> 0), for n> 1, where 

Z l n (t) = inf{s > 0\D n (s) + L n (s) > X n (t) + 5 n (t) + L n (t) - (1 - 1/n)} 
(3-6) _ 

= inf{s > 0\D n (s) + L n (s) > X n (0) + A n (t) - (1 - 1/n)}, 

Zl{t) = inf{s > 0|,D n ( S ) > X n {t) + D n (t) - (1 - 1/n)} 
(3-7) _ 

= inf{s > 0\D n (s) > X n (0) + A n {t) - L n {t) - (1 - 1/n)} 

for t > 0. For the next step, define the processes U l n = Z l n — e and U% = Z% — e 
for n > 1. Then we see that the bounds on our virtual waiting time process 

(3.4) and (3.5) can be written as V l n = {U l n ) + and V% = (U%) + , for n > 1. 
We would then want to use the corollary of [11], along with the assumptions 
that X = 1 and V = to get 

(3-8) c n ^ = c n (Zi-e) + ^£^, 

(3-9) c n V: = c n (Z--e) + ^^ F 

in D as n — > oo. However, notice that the right-hand side of the condition in 
(3.7) does not satisfy the conditions of the corollary. In particular, X n (0) + 
A n — L n — (1 — 1/n) is not necessarily a nondecr easing element of D for 
each n > 1. Therefore, we cannot deduce (3.9) directly. We can resolve this 
problem by using linear interpolations that bound the original processes 
from above. 

Notice from (3.7) that 

Zl = D- 1 o (X n (0) + A n -L n -(l- 1/n)) = b~ l o Y n , 

where we define Y n = (X n (0) + A n — L n — (1 — 1/n)), for n > 1. Since Y n 
is not necessarily a nondecr easing element of D, Z™ is not necessarily an 
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element of D (see Example 13.2.1 of [15]). Therefore, we construct linear 
interpolations -D" 1 and Y n of -D" 1 and Y n , respectively, in such a way that 

D- 1 ®^!)- 1 ® and Y n (t) < Y n (t) 

for all t > 0, and D~ l is nondecreasing, for each n > 1. The construction can 
be carried out as follows: when stepping down, linearly interpolate between 
the right endpoint of the previous step and the midpoint of the next step. 
Similarly, when stepping up, linearly interpolate between the midpoint of 
the previous step and the left endpoint of the next step. We show how to 
construct the interpolation in Figure 1. 

Then, for each n > 1, let Z% = D~ x o Y n , and notice that, since D^ 1 is 
nondecreasing, we must have Z™(t) < Z^(t) for t > so that 

(3.10) V:(t)<V:(t) fort>0, 

where = (Z%-e)+. 

By Lemma 7.1, the error caused by these linear interpolations is asymp- 
totically negligible. Also, by our assumed limit (3.2) and Theorem 13.7.2 of 
[15], 

^{D" 1 - D~\Y n -(A- L)) (- o D _1 ,A - in D 2 

as n-> oo, jointly with (3.1) and (3.2). Therefore, using the converging 
together lemma (Theorem 11.4.7 of [15]) gives us 

^(Dn 1 ~ D-\Y n -(A- L)) =► (- (B^ oD-\A-L) in D 2 

as n— > oo, jointly with (3.1) and (3.2). Since the limit process in the last 
limit belongs to C 2 , the limit holds in the stronger uniform topology. But 
since the prelimit processes also belong to C 2 the limit holds in C 2 . Using 
a version of the composition result Theorem 13.3.3 of [15] for elements of C 
with the continuous mapping theorem gives us 

c n V: = c n {Zl - e)+ = c n {D~ l oF„-e)+4£ in C 
as n — ► oo. Then, using the bounds (3.10) and (3.3), gives us our result. □ 

A similar bounding argument does not work in the ED case, because the 
limits of the bounding processes do not coincide. 

• e , «• 

• 1 © • ■- 



Fig. 1. Upper bound linear interpolation for step functions. 
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4. Limits from stopped arrival processes. In this section, we prove heavy- 
traffic limits for virtual waiting-times for G/G/n + G queues from associ- 
ated heavy-traffic limits for basic queueing processes with the arrival process 
stopped at some time r > 0. The results here will be applied to the Marko- 
vian M/M/n + M model in the ED limiting regime in Section 6.1. 

We present results for one-dimensional virtual waiting times and virtual- 
waiting-time processes, respectively, in the next two subsections. 



4.1. One- dimensional virtual waiting times. To get one-dimensional lim- 
its for the virtual waiting time at some fixed time r > 0, the idea is to stop 
the arrival process to the system at time r. We add the superscript r to 
all processes defined in Section 3 to denote the same processes but with the 
arrival process to the system stopped at time r. We then have for n > 1, 

VT{t) = inf{ S > 0\Dl(t + s)+ U n (t + s) > Xl(0) + A T n (t) - (n - 1)} 

for r,t > 0, and V n {t) = V*{t) for t > 0. 

For each r > 0, we assume we have the limits 

(4.1) (A T n ,D T n ,L T n ,X^^(A T ,D T ,L T ,X T ) in D 4 

as n — > oo, where A T , D T , 17 and X T are continuous deterministic functions 
and D T + 17 has continuous positive derivative {D T + L T )' . In addition, we 
assume that 

(4.2) (A T n ,Dl,L T n ,X^^(A T ,D T ,L T ,X T ) in D 4 

as n — > oo, where A T , D T , L T and X T are continuous processes. 

Then if we define for each r > the first-passage time processes = 
(j£(t),t>0) for n> 1, and Z T = {Z T (t), t > 0), where 

Zl{t) = inf{ S > 0|^(a) + Z;(8) > X;(0) + Alit) - (1 - 1/n)}, 
Z T (t) = inf{s > 0\D T (s) + Z T (s) > X T {0) + ^ r (t) - 1} 

for t > 0, the processes = Z^ — e for n > 1, and U T = Z T — e, then we 
see that for each r > and n > 1 our virtual-waiting-time processes can be 
written as V£ = (0^) + . We first prove a limit for the processes and then 
use this limit to prove our virtual- waiting-time heavy-traffic limits. For each 
r > 0, define the process U T = {U T (t),t > 0) by 

(4.3) fjrf t) = ^n(O) + i T (t) - £> T (^ r (t)) - ^ T (^ r (i)) 



(,D T + L T ) / (Z r (t)) 



Note that [/ r is a continuous process. 



20 R. TALREJA AND W. WHITT 

Lemma 4.1. For each r > 0, under the assumptions (4-1) and (4-2), 
c n {Ul - tJ T ) ^U T inD 
as n — > oo, jointly with the limits (4-1) and (4-2), where U is given in (4-3). 

Proof. Applying the corollary of [11], we get 

Cn{Ul - U T ) = c n [{Zl - e) - (Z T - e)\ = c Tl {Z T n -Z T )^U 
in D as n — ► oo. □ 

Theorem 4.1. For each r > 0, under the assumptions (4-1) and (4-2), 
we have: 

1. IfX(r) = 1, then 

K(r) = c n U-(r)+ U t (t)+ = (f *}% ( m R 

{D + Ly{r) 

as n— > oo, jointly with the limits (4-1) and (4-2). 

2. If X{t) > 1, then 



V n {r) = c n (^(r)+ - U t (t)) U^t) 



in 



as n— >oo, jointly with the limits (4-1) and (4-2), where U t {t) is given 
in (4.3). 

Proof. Part 1 follows immediately from Lemma 4.1 and the fact that 
Z t {t) = r, X t (t) = X(r), D t (t) = D{t) and Z r (r) = L(r) for all r > 0. 

For part 2, notice that the condition X{r) > 1 implies U t (t) > 0. By 
Lemma 4.1, we get the limit 

UI{t) =► U t {t) in R 

as n — > oo. Therefore, for each e > and < 5 < U t (t), we have 

n\cn{tJ T n {r) + - U t (t)) - CnitJKr) - U T (r))\ < e] 
= P[ Cn |f7 n >)+-C7;( T )|< e ] 
(4.4) >PK(r)>0] 

>F[(%(t)>Ut(t)-6\ 
>F[\U^t)-U-(t)\<5]^1 

as n — > oo. Combining this limit with the convergence in Lemma 4.1 and 
Theorem 4.1 of [2], we complete the proof. □ 
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4.2. Virtual waiting time process. To get an associated limit for the 
virtual-waiting-time processes, we again assume that the limits (4.1) and 
(4.2) hold, but where the queueing processes are understood to be processes 
in Djj indexed by both r and t. To be precise, for each n > 1 we define the 
processes 

= ( X (2) (T) t ) > Q) ^ A (2) _ ( A (2) ( Tj t ) s r> t > 0) , 

Z?i 2 ) = (2#> (r, t ) , r, t > 0) , L^ = (L^(r,t),r,t>0), 

so that 

X«(T,t)=^(t), 4 2 )(r,t)=A;(t), 
D®(T,t) = I%<t), L^(r,t) = L T n (t) 
for r, i > and n > 1. Next, define the scaled processes 

V(2) 4 (2) n (2) r (2) 

-^(2) = £n_ J(2) = flL_ £){2) = t^_ TW = — 

n ~ n ' n ~ n ' n ~ n ' n ~ n ' 

and 

*( 2 ) ee Cn (X( 2 ) - X( 2 )), i( 2 ) ee c„(4 2) " ^ (2) ), 

££)^(£>£) -£>«), 4 2) -c n (Zf-Z( 2 )), 

where c n — > oo and n/c n — > oo as n — > oo, andX( 2 ), A< 2 >, L»( 2 ) and are 
given by 



(4.5) 



xW(r,t) = X T (t), A^{r,t) = A T {t), 



Note that , £>( 2 ), L( 2 ) and X^ 2 ) are continuous deterministic functions 
and, for each r > 0, Z)( 2 )(t, •) + L^ 2 \t, •) has continuous positive derivative. 
Then we assume we have the limits 

(4.6) (AP) t D&>,L® t XV>) => (^ 2 ),I)( 2 ),Z( 2 ),X( 2 )) in 
as n — ► oo, and 

(4.7) => (i (2 \^ (2 U (2 U (2) ) in 0* 
as n — > oo, where X( 2 \ i^ 2 ), L>( 2 ) and Z( 2 ) are given by 

x( 2 \r,t) = X T (t), A^(r,t)^A T (t), 

(4.8) 

D®( T ,t) = D T (t), m{r,t) = L T {t). 
Note that X( 2 ),i( 2 ),Z»( 2 ),Z( 2 ) G C c . 
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— (2) — (2) 

Next, define the processes U n = (U n (t, t),r, t > 0), for n > 1, and the 
process C> (2) = (C> (2) (r, i), r, t > 0) so that 

(4.9) U,P(r,t)^U^t), uW(T,t) = UT(t) 
for r, i > and n > 1. Then, as in Lemma 4.1, we have 

Lemma 4.2. Under the assumptions (4-6) and (4- r 0, we have 

(4.10) c n (^ 2) -i/ (2) )^?7 (2) iiiDo 

as n — ► oo jointly with the limits (4-6) and (4-1), where and are 
given by (4-9). 

Proof. The proof follows the proof of Lemma 4.1, except instead of 
using Corollary 1 of [11] we use Theorem 2.9 here. □ 

Now define the projections U° = {U°{t),t> 0), for n > 1, and U° = (U°(t),t> 
0) so that 

u°(t) = uP(t,t), (%(t) = uW(t,t) 

for t > 0. Then we have 

Theorem 4.2. Under the assumptions (4-6), (4-V and X^ eCc, we 
have the following: 

1. // X(t) = 1 for all t > 0, then 

Vn = C n (U° n ) + (U°) + = jJ^y in D 

as n— > oo, jointly with the limits (4-1) and (4-%)- 

2. // inf t > X(t) > 1, then 

Vn = c n ((U° n ) + - U°) ^U° inD 
as n — > oo, jointly with the limits (4-6) and (4-1)- 

Proof. The limit in Part 1 follows immediately from the limit (4.10) by 
an application of the continuous mapping theorem with the projection map 
of Lemma 2.3. The second equality then follows from the fact that Z l (^) = t, 
X\t) = X(t), D\t) = D(t) and Z*(t) = L(t) for all t > 0. 

For Part 2, notice the condition inf r >oA'(r) > 1 implies inf T >o U t (t) > 
so that as in (4.4) we get for all e > and T > 0, 

n\K((Un) + - U°) - c n (U° - U°)\\ T < e] 1 as n - oo. 

Combining this limit with the limit (4.10) and the continuous mapping the- 
orem with the projection map of Lemma 2.3 completes the proof. □ 
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5. QED. In this section and the next, we apply the results of the pre- 
vious two sections to prove heavy-traffic limits for virtual waiting times in 
the Markovian M/M/n + M model in the QED and ED limiting regimes, 
respectively. We assume that inter arrival times, service times and patience 
times after time are mutually independent and independent of the initial 
conditions. To construct our heavy-traffic limits, we consider a sequence of 
M/M/n + M models indexed by n > 1. The nth model has n servers, each 
with service rate p, arrival rate A n , and individual abandonment rate 9. The 
QED limiting regime of [5] is characterized by the limit 

lim v/n(l — p n ) = 0, — oo < 8 < oo, 

n— >oo 

where p n = X n /np is the traffic intensity in the nth model. Let the basic 
processes associated with the models, i.e., X n , A n , D n , L n , V n , X n , A n , D n , 
L n , V n , X n , A n , D n , L n and V n , for n > 1, be defined as in Section 3. For 
this QED case, under initial conditions such that X n (0) 1 in K as n — > oo, 
we will have the fluid limits 

(5.1) X = 7], A = ue, D = pe, L = 0, V = 0. 

Since X = n and L' = 0, the results of Section 3 apply. 

We start with the following theorem (Theorem 2 of [3]) for all processes 
except for V n . See Theorem 7.1 of [9] for an alternative proof. 

Theorem 5.1. Consider the sequence of M/M/n + M models defined 
above. If X n (Q) X(0) in R as n — > oo, then we have the fluid limit 

(5.2) (A n ,D n ,L n ,X n )^(A,D,L,X) inD A 

as n^ oo, where the limit is given in (5.1). Furthermore, if X n (0) X(0) 
inM as n — > oo, then we have the diffusion limit 

(A n ,D n ,L n ,X n )^(A,D,L,X) mD 4 

as n — > oo, where X is the unique solution to the stochastic differential equa- 
tion 

dX(t) = -fi/3 - p(X(t) A0)dt- 9{X{t) V0)dt + y/JidBi(t) - ^/J2dB 2 (t), 
and 

A{t) = yfjiB X {t), 

D(t) = ^J2B 2 (t) + p f\x(s) A 0) ds, 
Jo 

L(t) = 9 f\x(s)V0)ds 
Jo 

for t>0, where (Bi,B 2 ) is standard Brownian motion in two dimensions 
and is independent of X(0). 
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Then, applying Theorem 3.1, we obtain: 

Theorem 5.2. Under the assumptions of Theorem 5.1, 

X+ 

V n in D 

as n — > oo, jointly with the limits in Theorem 5.1. 

We can use Theorem 5.2 to get a heavy-traffic limit for potential waiting 
times of customers in the queue (see [13]). Let W n = (y/nWn^ +1 ,t > 0), 
where denotes the potential waiting time of the ith customer to enter 
system n, for i > 1, n > 1. Then using the fact that 

W \nt\+l = y n{ r r \nt\+l^ t > Q) 

where is the arrival time of the ith customer in the nth system, for i > 1, 
n > 1, along with the continuous mapping theorem with the composition 
map, we get 

Corollary 5.1. Under the assumptions of Theorem 5.1, 

as n— > oo, jointly with the limits in Theorem 5.1. 

Finally, we prove convergence of steady-state virtual waiting times. The 
argument is given in the proof of Theorem 3 of [3] , but we repeat it here for 
completeness. 

Corollary 5.2. There exist random variables X(oo) and V n {oo) such 
that 

(5.3) X(t)=>X(po) in R as i->oo, 

(5.4) V n (t) =>V n (oo) in R as t — > oo 

for n > 1. Furthermore, under the assumptions of Theorem 5.1, 

- . , X(oo) + 

V n (oo) => — — — in R 

as n — > oo. 

Proof. We first recall from [3] that X{oo) is indeed well defined. Also, 
recall from Theorem 2* of [3] that we can interchange the heavy-traffic and 
steady-state limits of X n (t). We can use this fact to deduce the interchange 
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for V n {t) as follows. Note that each of the process X n , n > 1, possesses 
a stationary distribution since it is simply a birth-death process that can 
easily be seen to be positive recurrent. For each n > 1, define X n (0) to have 
this distribution so that we have a stationary version of the process X n . 
Because of the birth-death structure of the processes X n , n > 1, for each 
t > and n > 1, V n (t) can be represented in terms of X n (t) by 

(X„(t)-n)+ 

V n (t) = X i,n, 
i=l 

where n represents the time between the (i — l)th and ith departure 
(either an abandonment or a service completion) from the queue in model 
n > 1. Therefore, for each n> 1, V n must also be a stationary process so 
that V n (t) has the same distribution for all < t < oo. Thus, for each n > 1, 
the random variable V^(oo) defined in (5.4) is well defined with the same 
distribution as V n (t), for any fixed t > 0. Projecting the result of Theorem 
5.2 onto a fixed t > and using the fact that we can interchange the heavy- 
traffic and steady-state limits of X n (t) then gives us our result. □ 



6. ED. We now study waiting-time asymptotics in the ED limiting regime 
using the results of Section 4. We again have a sequence of M/M/n + M 
models with n servers in model n > 1 , but now customers arrive to the sys- 
tem at the rate A n = nX, for some A > 0. In all the models, each server has 
service rate (i < A and each customer has abandonment rate 9. As in Section 
5, we assume that interarrival times, service times and patience times after 
time are mutually independent and independent of the initial conditions. 
Again, let the basic processes associated with the models, i.e., X n , A n , D n , 
L n , V n , X n , A n , D n , L n , V n , X n , A n , D n , L n and_V n , for n > 1, be defined 
as in Section 3. Under initial conditions such that X n (0) =4> 1 in R as n —> oo, 
we will have the fluid limits 



(6.1) 



X = (q+i)r), q=^Y~~> A = Xe, D = fie, 
L = 0qe = (A — u)e, V = wri, w = -ln(— 

e \fi 



for t > 0. 

We start by stating ED fluid and diffusion limits for A n , D n , L n and X n 
These appear in Theorem 2.1 and Corollary 2.1 of [17]. 



Theorem 6.1. Consider the sequence of M/M/n + M models defined 
above. If X n (0) =>■ X(0) in R as n — > oo, then we have the fluid limit 

(6.2) (A n ,D n ,L n ,X n ) (A,D,L,X) in D 4 
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as n — > oo, where the limit is given in (6.1). Furthermore, if X n (0) =^ A"(0) 
in M. as oo, then we have the diffusion limit 

(6.3) (A n ,D n ,L n ,X n )^(A,D,L,X) in D 4 

as n — > oo, where X is the unique solution to the stochastic differential equa- 
tion 

dX(t) = -9X(t)dt + V\dB 1 (t)- ^JIdB 2 (t) - ^/\^JidB 3 (t), 

and 

A^^VXB^t), 
t)(t) = v^B 2 (t), 

L(t) = y/X^B^t) + e f X(s) ds 

Jo 

for t>0, where (Bi, B2, B 3 ) is standard Brownian motion in three dimen- 
sions and is independent of X(0). 



6.1. Diffusion limit for virtual waiting time. We now prove a diffusion 
limit for the virtual waiting-time process in the ED regime, Theorem 6.3. 
As discussed in the Introduction, our result is a special case of Theorem 
4.1 of [7], but they do not provide a complete proof. We first state the 
following result, Proposition 6.1, which gives us the limits (4.1) and (4.2) 
in the case of the M/M/n + M model in the ED limiting regime. It is a 
special case of Theorems 3.1 and 3.2 of [7], which are proven using strong 
approximation (see also [6]). The result is also a special case of Theorems 2.1 
and 2.2 of [12], which are proved using a different approach. We will then 
combine Proposition 6.1 with Theorem 4.1 to obtain the one-dimensional 
heavy-traffic limits for virtual waiting times, Theorem 6.2 below. Afterward, 
we obtain heavy-traffic limits for virtual waiting time processes in Theorem 
6.3. 



Proposition 6.1. Consider the sequence of M/M/n + M models de- 
fined above and fix r > 0. // X n (0) X(0) in M as n — > 00, then we have 
the fluid limit 

(A T n ,Dl,L T n ,Xl) (A T ,D T ,L T ,X T ) m D 4 
as n — > 00, jointly with the limit (6.2), where 
Xe -e{t-r)+ _ ^ 



(6.4) X7(t) 



+ 1, t<T + W, 

- M (t-( T +tt)) t > T + w 
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(6.5) A T (t) = A(tAr), 

(6.6) D T (t) = n(t A (r + w)) + — (1 — e -/*(*-(T+«)) + ) } 

A* 

(6.7) U{t) = (A - M )(t At) + ^(1 - e^^o ) - M [t - T ]« 

/or i > 0, where w = (1/0) ln(A//i). Furthermore, if X n (0) => X(0) in R as 
n — > oo, i/ien we /iai>e i/ie diffusion limit 

{A T n ,b r n ,Ll,X T n ) (A\D\L T ,X T ) m D 4 

as n —> oo, jointly with the limit (6.3), where X T is the unique solution to 
the stochastic integral equation 

X T (t) = X T (0) - j\fil {s > T+w} X T (s) - 81 {s<T+w} X T (s)} ds 

(6.8) + V\Bx(t At) - ^J2B 2 (j\x T (s) A 1) ds^j 

-Ve~B 3 n\x T {s)-i) + ds^, 

and 

(6.9) A T {t) = \f\B 1 {t At), 

(6.10) D T (t) = nj^ l {s > T+w} X T (s) ds + ^J1B 2 (J q {X t (s) A 1) ds^j , 

(6.11) L T (t) = j* l {s<T+w} X T (s) ds + V9B 3 (jj^ W - 1) + ds) 

for t>0, where (B\, B 2 , B%) is standard Brownian motion in three dimen- 
sions and is independent of X(0). 

We omit the proof of Proposition 6.1 because it is a direct consequence of 
the two-parameter version in Proposition 6.2 below. The distribution of X T 
in (6.8) is complicated because of the deterministic time transformation in 
the Brownian motions, but existence and uniqueness follows from Theorem 
2.11. Combining Proposition 6.1 with part 2 of Theorem 4.1, we immediately 
obtain the following one-dimensional result. 

Theorem 6.2. Under the assumptions of Theorem 6.1, for each r > 0, 
as n— > oo, jointly with the limits in Proposition 6.1. 
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We now establish a two-parameter version of Proposition 6.1. 

Proposition 6.2. Consider the sequence of M/M/n + M models de- 
fined above. If X n (0) =^ X(0) in R as n — > oo, then we have the fluid limit 

(AW , W , Z( 2 ) , ) (A^ ,DW, IW ,XW) in D% 

as n — > oo, jointly with the limit (6.2), where X^i ; A( 2 \ D^) and LW 
are given by (4-5) and (6.4)-(6.7). Furthermore, if X n (0) X(0) in R as 
n — > oo , then we have the diffusion limit 

(AW,f)W,LW,xW) -> (AW,dW,lW, X W) ^n D% 

as n — > oo, jointly with the limit (6.3), where X^ , A^ , E>W and iP^ are 
Sivera 6y (4.8) and (6.8)-(6.10). Furthermore, (A^\D^\W,X^) E C£. 

Proof. The proof follows the same steps as the proof of Theorem 7.1 
of [9] except some extra care must be taken to apply the continuous map- 
ping theorem with our integral mapping. First, after centering, the system 
equation becomes 

Af ) (t, t) = XW ( T) 0) + Mi 2 } (r,t)- M$ (r, t) - Af« (r, t) 

(6.12) +\ n (tAr)- n (XP(r,s) An)ds 

Jo 

-6 [ t (xW(r,s)-n) + ds 
Jo 

(2) (2) 

for r,t > 0, where the processes M„ ^ = (M^{ (r, t),r,t > 0) are defined for 
* = 1,2,3 by 

M<3(r, i) = A(\ n (r A i)) — A n (r A t), 

M® (r, t) = S (/x £ (X® (r, s) A n) da) - /x jf ' ( X^ (r, a) A n) da, 

(r, t) = L (0 J* (XW ( T , s ) - n ) + ds^j - 9 j* (X^ (r, a) - n) + ds 

for r, i > 0, and A, S and L are the unit Poisson processes generating arrivals, 
departures and abandonments in the setting of Theorem 6.1. 

(2) 

Fluid limit. Since for each n > 1 we have Xn (r, a) < X n (s) for all r, a > 
0, we have as in Lemma 4.5 of [9] 

M (2) 

(6.13) _^L^o in D D w.p. 1 fori = 1,2, 3, 
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jointly as n — ► oo using the FSLLN for Poisson processes. We only apply the 
weaker FWLLN consequence. 

Then, using the integral map of Theorem 2.12 with the function h\ : 
E x [0, oo) 2 -> R defined by 

hi(x,s,t) = -l{ t > s+w }Hx - l{ t<s+w }0(x - 1), x G R,s,t > 0, 

X T must satisfy the integral equation 

X {2) (r, t) = X (2) (r, 0) + A(r At) - //((r + w) A t) 

rt _ 
- J Q l{ s >r+«,}^ (2) (r, a) + l {s<T+ „, } #(X (2) (r, s) - 1) ds, 

for r, t > 0. Lemma 7.2 shows that h\ satisfies the conditions of Theorem 
2.12. It is easy to see by inspection that this integral equation has solution 

— (2) - (2) 

given by (6.4). We can now use (6.13) to get the limits for A n , D n and 
L n jointly. 

Diffusion limit. Using our two-parameter integral mapping (Theorem 
2.10) and two-parameter composition mapping (Theorem 2.4) with the con- 
tinuous mapping theorem, we get the limits 

M (2) 

(6.14) M n 2 ] = — ^ in D D for i = 1,2, 3, 

\/n 



jointly as n — > oo, where 

^ 2) (r,f) = Bi((A(rAt))), 

M 2 (2) (r, t) =B 2 (fJiJ (X {2) (r, s) A 1) ds 
M 3 (2) (r, t) = B 3 (ej o (X^ (r, S ) - 1) + ds) . 

Notice that 6 Cc so that the applications of the continuous mapping 
theorem with Theorems 2.10 and 2.4 are justified. 

By Lemma 7.4 and system equation (6.12), there exists a sequence of 
processes H n = {H n {r, t),r,t>0), n> 1, such that 

(6.15) H n ^0 in Do as n — > oo 

and 

Af ) (r, t ) = Af) (r, 0) + M$ (r, t) - M$ (r, t) - Afg (r, t) 
+ A n (r At) - a*((t + w) At)n + ^H n (r,t) 

(6.16) 
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- fJ> J l{s>T+w}X ( nXT,s)<ls 
rt 

-0 J 1{s<t+w}(xP(t,s) -n)ds. 

Centering X n 2 \r, t) by nX^ (r,t) in (6.16) and dividing by y/n, we get 
XP ( T , t) = X^ (r, 0) + Mi 2 } (t, t) - J&g (r, t) - (r, t) 

(6.17) +H n (T,t)-fi l {s>T+w} X n 2 \ T ,s)ds 

Jo 

ft 

~° Jo 1 { S<T + U, } X " ( T ' S ) ds 

for i > 0. 

Define the function ho : R x [0, oo) 2 — > R by 

/i (x, s, i) = -l {t > s+tu }pia; - l{ t<s+tu }6>x, x G R, s > 0. 

It follows from Lemma 7.2 with a = that this function satisfies the condi- 
tions of Theorem 2.12. Applying the continuous mapping theorem with the 
integral map of Theorem 2.12 and the addition map of Theorem 2.2 to the 
limits (6.14) and (6.15) gives us our limit for X^. The limits for A n and 

* (2) * (2) * (2) 

L n ' follow directly from the limits for M n [ and M n g . Then, the limit for 

D n follows from (6.17) and the continuous mapping theorem with addition. 

Theorems 2.4, 2.10 and 2.12 also show that the property of being in Cc 
is preserved with respect to applications of these two-parameter mappings. 
□ 

Now define the process X w = (X w (t),t> 0) so that 

X w {t)=X (2 \t,t + w). 

Then, combining Proposition 6.2 with part 2 of Theorem 4.2 and simplifying 
we obtain the following theorem. 

Theorem 6.3. Under the assumptions of Theorem 6.1, we have 

X w 

V n => in D 

A* 

as n — > oo, jointly with the limits in Theorem 6.1. 



Finally, paralleling Corollary 5.1, we get the following corollary. 
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Corollary 6.1. Under the assumptions of Theorem 6.1, 
cis n — > oo, jointly with the limits in Theorem 6.1. 



6.2. Steady-state virtual waiting time. We now prove convergence of steady- 
state virtual waiting times under diffusion scaling in the ED limiting regime. 
It is interesting that, although we are only proving a one-dimensional limit, 
our proof uses convergence in D. 

Just as in Corollary 5.2, we can define random variables V^(oo), for n > 1, 
in the ED setting. From (2.26) of [17], we know that 

(6.18) V n (oo) ^w = ^ ln Q) inM 

as n — > oo. This result is obtained there by analyzing an ODE for the fluid 
limit of the queue-length process with arrivals turned off. We now give a 
different argument for (6.18) as well as a refinement. Let Q n {oo) denote the 
steady-state queue-length (excluding customers in service) in system n > 1. 
As in the proof of Corollary 5.2, because of the birth-death structure of Q n , 
V n (oo) can be represented directly in terms of Q n (oo) by 

Qn(oo) 

(6.19) V n (oo)= X i,n, 

i=l 

where JQ„ represents the time between the (i — l)th and ith departure 
(either an abandonment or a service completion) from the queue in model 
n > 1 . For each n ^ 1 , the Xi n are independent of each other and of Q n (co), 
and each Xi >n is exponentially distributed with mean l/(nfi + i9). 
First, we prove the following lemma. 



Lemma 6.1. Let the triangular array (^Q, n ) be defined as in (6.19). Then 

/\n-\ \ 

\fn Xi jV , — c \ ^ B o d in D 



as n — > oo , where 



1 / 9t N 
c(t) = -]n 1 + -), d(t) 



[i ) ' + Ot) 

for all t > 0, and B is standard Brownian motion. 
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Proof. It can easily be verified that the asymptotic mean and variance 
of the sum in (6.19) are c(i) and d(t)/n, respectively. Also, the triangu- 
lar array (Xi tTl ) satisfies Lyapunov's condition. Therefore, the functional 
Lindeberg-Feller central limit theorem applies to the partial sum sequence 
and we get our result; see Theorem 4.1 of [10]. □ 

We now give our main result of the section. Note that one could erro- 
neously conclude from Theorem 6.3 that the limiting distribution of V^(oo) 
has variance (l/// 2 ) Var(X(oo)) = \/{]i 2 6). We see here that this is in fact 
not the case. What is happening is that the numerator of the limit at a fixed 
t > 0, (t, t + w), corresponds to the diffusion-scaled number-in-system 
observed w units of time after stopping the arrival process. This delay in 
observing the number-in-system is crucial to the analysis and is the reason 
one can not simply use the steady-state distribution of X. 

Theorem 6.4. Under the assumptions of Theorem 6.1, 



as n — ► oo, where V^(oo) = yfn(y n {po) — w). 

Proof. Let Q n = (Q n (t),t > 0), where Q n (t) is the number of customers 
in the queue (excluding customers in service) at time t in model n>l. Note 
that for each n, the process Q n possesses a stationary distribution since 
Q n = (X n — n) + and X n is simply a positive recurrent birth-death process. 
Redefine Q n (0) to have this distribution so that we have a stationary version 
of the process Q n . By Theorem 2.3 of [17], we have Q n (0) =4> q in E as n — ► oo, 
so that applying Theorem 6.1 and the continuous mapping theorem gives 
us the fluid limit for the queue length process Q n (t) qr], in D as n — ► oo. 
Joining this convergence with the convergence in Lemma 6.1 using Theorem 
11.4.5 of [15] gives us 



as n — > oo. Projecting the second component on to a fixed t > gives us 




M 




in D 



,2 



\ \t=0 / / 

Applying Proposition 13.2.1 of [15] gives us 




as n 



oo 



(6.20) 




m 



M 
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as n — > 00. Then, 

/Q n (t) \ /Qn(t) \ 

E - c («)J = ^ ( E x *>- - c (0n(t))J 

+ V^(c(Qn(i))-c(g)) ) 

where c(g) = u>, for w in (6.18). To compute the limit of the second term 
above, we use the Taylor approximation 

c(Qn(t)) = c(q) + c'(q)(Q n (t) -q)+ c"(q)(Q n (t) - qf + • • • 

= c(q) + c'(q)(Q n (t) -q)+ o ¥ (n~ 1 ^), 

where < e < 1/2. The second equality above follows from the continuous 
mapping theorem and the fact that 

n 1/2+e (Q n (t) - qf = n- x l**{yfc{Q n {t) - q)f ^0 in R 

as n — > 00. Observing that c'(q) = 1/A gives us 

(6.21) x/E(c(Q n (t))-c{q)) = c'(q)Q n (t) + o ]? (n- e )^jQ(t) in E 

as n — > 00. Combining the limits (6.20) and (6.21) and using the continuous 
mapping theorem with addition we get 



_/Qn(t) \ -. 

V 7 ^ E *v» ~ c(q)j => B(d(q)) + -Q{t) 



m 



as n — > 00. Finally, since V^(oo) is distributed as the stationary distribution 
of the left-hand side above, and the stationary distribution of the limit is 
normal with mean and variance 

%) + ;pVa r( Q) = A + sr_L, 

we get our result. □ 

7. Additional lemmas. We use the following lemma to show that the 
error caused by our linear interpolations in the proof of Theorem 3.1 is 
asymptotically negligible. 

Lemma 7.1. For the linear interpolations D^ 1 and Y n of D~ l and Y n , 
respectively, we have for all T > 0, 



V^\\(D-\Y n )-(D-\Y n )\\ T ^0 



in 



as n — > 00, jointly with the limits in Theorem 5.1. 
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Proof. Notice that for each T > our construction gives us 

(7.1) \\D- 1 -D- 1 \\ T = J T {D- X ) and \\Y n -Y\\ T = J T (Y n ), 

which both converge to as n — ► oo since D~ l and Y n are continuous by 

(5.2) and Corollary 13.7.3 of [15]. Therefore, by (7.1), we have 

sft\\{b-\Y n )-{D-\Y n )\\ T 

n sup \\(D-\t)-D-\t),Y n (t)-Yn(t))\\oo 

0<t<T 

n sup maxUD-^t) - D- 1 ^,^) -Y n {t)\} 

0<t<T 

= v / nmax{||£'~ 1 - D^Wt, \\Y n - Y n \\ T } 
= ^max{J T (D- 1 ),J T (F ri )} 

= max{ J T {V^(D- 1 - D- 1 )), J T (^(Y n -{A- £)))} => 

in R as n — > oo, since the processes D and A — L are continuous. Here, we 
assume the maximum norm || • \\oo on R without loss of generality. D 

The following lemma shows that the regulator map used in Proposition 
6.1 satisfies the conditions of Theorem 2.12. 

Lemma 7.2. For o£l, the function h a :Rx [0, oo) 2 — > R defined by 
h a (x, s, t) = -l {t > s+w} fix - l {t<s+w} 6(x -a), x G R, s, t > 0, 
satisfies the conditions of Theorem 2.12. 

Proof. h a clearly satisfies condition 1 of Theorem 2.12. We now show 
h a satisfies condition 2 of Theorem 2.12. Fix T > 0, s\,s 2 > 0, xi,x 2 € -Dd, 
and an increasing homeomorphism A on [0, T] having strictly positive deriva- 
tive A'. We then have 

t 

h a (xi(s 1 ,u),s 1 ,u) - h a (x 2 (s2,X(u)),s 2 ,X(u))\ du 
t 

|-l{ u > Sl+lu }/xxi(si,n) - l{ n<S2+w }6 , (xi(si,n) - a) 

+ ^{X(u)>s 2 +w}^2(s 2 , X(u)) + I{\( u )<s 2 +W}Q(X2(S2, X(u)) - o)| du 

-vf ^{u^si+wyxiisuu) - l {x( , u) > S2+w} x 2 (s2,X(u))\du 
+ 8 [ l 1 {u<«i+w} x l( s l>' li ) _ 1 {\{u)<s 2 +w}X2(s 2 ,X(u))\du 
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+ Jo ^t u<Sl+w } ~ ^{Hu)<s 2 +w}\ du 



< 



[1 I \l{ u > si +w} -^-{X(u)>S2+v,}W x 2(s2,X(u))\du 



+ V / l{u> Sl +w}\xi(si,u) - x 2 (s 2 ,X{u))\du 



+ 6 ^{uo-i+w} - t \\(u)<s 2 +w}\\x2(s2, A(n))| du 



+ 9 / l{ u < Sl +w}\xi(si,u) - x 2 (s2,X(u))\du 
Jo 

Jl ~^ a J Q \l{u< Sl +w} - ^{\(u)<s 2 +w}\du 

< li\\x s 2 2 \\ t + ^(||^ 2 || T + a)|A- 1 ((a 2 + w) A T) - (( Sl + w) A T)\ 

+ (n + 9) \xi(sx,u) - x 2 (s 2 , A(u))| du 
Jo 

<fi\\x s 2 2 \\ T + 6{\\x s 2 2 \\ T + a) 

x (\X~\(8 2 + w) A T) - \- X {( Sl + w) A T)\ + HA" 1 - e\\ T ) 

+ (fj, + 9) \xi(sx,u) — x 2 (s 2 , A(u))| du 
Jo 

< mII^IIt + 0(\\x?\\ T + o)||(A- 1 )'|| r (|a2 - «i| + ||A - e|| T ) 

+ + 6>) / |zi(si,u) -a; 2 (s 2 ,A(u))|du, 
Jo 

where the last inequality follows from the mean value theorem. Since 

iKA-'yiir < : - r 1 VM <°°, 

info< s <TA'(s) 

this completes the verification of condition 2 of Theorem 2.12. □ 

We now prove the approximations used in the proof of Proposition 6.2. 
First, we have the following lemma, which is essentially a two-parameter 
version of Theorem 10.1 of [6]. 

Lemma 7.3. In the context of Section 6, if X n (0) =>• X(0) in M as n — > oo, 
then for all S, T > 0, 



||X< 2) -.? <2, |I S .T = 0P(-^ 



1/2-e y " 
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Proof. For each n> 1, define the process Y n = (Y n (T,t),T,t > 0) by 

A n , 



Y n (r,t) = \X n (0)-X(0)\ + 



n 



-{t At) - X(t At) 



+ -|A(A n (tAT))-A n (tAT)| 



+ - 
1 

+ - 



S^riu J\x^(t, s) A 1) ds^ - nn J\x^\t, s) A 1) ds 



R\nn j (xP(T,s)-l) + ds 
-nn \ (XW(t,s) - l) + ds 



Then we have for all e > 

(7-2) \\Y n \\s,T 



n 



1/2-e l> 



by the assumed limit X n (0) =4> X(0) in IR as n — > oo, the assumption A^ — 
nX, the limits (6.14), and Theorem 2.2. Now by the integral representation 
(6.12), for each n > 1 we have 

\xP(r,t)-X^(r,t)\ 



H(XW (t, s) A 1) + #(X n 2 > (t, S ) - 1)+ 
/i(X( 2 ) (t, s) A 1) + d{X^ ( Tj s) - 1)+ ds 



<Y n (r,t) + 



<Y n (T,t) + (u + 6) / |X( 2 )(T, S )-Z( 2 )(T, S )|d S 



for all t, t > 0, where the second inequality follows from the fact that the 
function 

h( s ) = -u(sAl)-9{s-l) + 
is Lipschitz. This implies 

II*? " * (2) lls.T < \\Yn\\s,T + {n + e) - X^ ^ ds. 

Jo 

Applying Gronwall's inequality then gives us 

ll*n 2) -* {2) H5,T<^ +e)T rn||5,T. 

The result now follows from (7.2). □ 



We use the following approximation in (6.16). 
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Lemma 7.4. Under the assumptions of Proposition 6.1, for all S,T > 0, 



sup 

0<T<S 

0<t<T 



H / (X^(t,s) An)ds + 9 / (xP(t,s) - n) + ds 



nt 

(7.3) - /z((r + w) A t)n - jf l {s > r+u , } /iX( 2) (r, s) 



- 1 {s<r+«)}^(^ 2 ' ) (t, s) - n) <is 



op(\/n). 



Proof. The key is to show that for each r > and e > the first and 
last times Xn {t, •) hits n are within op(l/n 1 / 2 ~ e ) of r + u> (see Figure 2). 
Fix T > t + w. By Lemma 7.3, for all 61,62 > 0, there exists N > such that 
for all n > N 

(7.4) PL4£]>l-<5 2 , 
where is the event 

On A 5 ^ we have for all < r < S, 

(7.5) inf{ S >0|X( 2 )(r, S )<n}> J Br, 



n(l + 5) 




T T + W t 

Fig. 2. The process deviating from the process nX T . 



38 

(7.6) 
where 

(7.7) 
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sup{s>0\XP(T,s)>n}<B T u , 



BJ = mi\ s > 



Xe -8(s-ry 



9 



+ l)_ 5in l/2+e < 



n 



(7.8) B T U = inf{s > 0\ne~^ s ^ T+w ^ + + (W /2+e < n}. 

It can then be easily checked by manipulating the inequalities in (7.7) and 
(7.8) that 



sup \BJ — (r + w) 

0<T<S 

SU P \ B l-( T + w ) 

0<T<S 



1 



n l/2-eJ> 
1 



n 



1/2-e 



We now show that 

rt 



sup 

0<T<S 

0<t<T 



fi / (xW(T,s)An)ds 



o 



/•t 

/x((r + w) M)n-n J l {s > T+w} X^ (t, s) ds 



The approximation for the second integral in (7.3) can be shown using the 
same procedure and this will complete the proof. Let < e < 1/2. By (7.5) 
and (7.6), for all c > there exists N' large enough so that for all n > N' , 
(7.4) holds and on the event Af^ we have for all r, t > 

f (X^ (r, s) A n) ds - (((t + w) A t)n + J l {s > T+w} X^ (r, s) ds^j 
f (X( 2) (t,s) An) (is 

-4 2 )(r, S )) + & 

(X^(r,s)-n) + ds 
\X^(r, S )-nX^s)\d S 

J T+w—c/n 1 / 2- ^ 

(■r+w+c/n 1 / 2 -* 
+ \XP(T,s)-nX T (s)\ds 







< 





1 { S <T+i«} ndS + 

+io— c/n 1 / 2- e 



T+W 
T+W 
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n l/2+e 

<c5i— jr— <c5 in 1/2 . 

By (7.4), the above holds with probability greater than 1 — fe, giving us our 
result. □ 
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